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We analyze numerically the magnetorotational instability of a Taylor-Couette flow in a helical 
magnetic field [helical magnetorotational instability (HMRI)] using the inductionless approximation 
defined by a zero magnetic Prandtl number (Pm — 0) . The Chebyshev collocation method is used to 
calculate the eigenvalue spectrum for small amplitude perturbations. First, we carry out a detailed 
conventional linear stability analysis with respect to perturbations in the form of Fourier modes 
that corresponds to the convective instability which is not in general self-sustained. The helical 
magnetic field is found to extend the instability to a relatively narrow range beyond its purely 
hydrodynamic limit defined by the Rayleigh line. There is not only a lower critical threshold at 
which HMRI appears but also an upper one at which it disappears again. The latter distinguishes 
the HMRI from a magnetically modified Taylor vortex fiow. Second, we find an absolute instability 
threshold as well. In the hydrodynamically unstable regime before the Rayleigh line, the threshold 
of absolute instability is just slightly above the convective one although the critical wavelength of 
the former is noticeably shorter than that of the latter. Beyond the Rayleigh line the lower threshold 
of absolute instability rises significantly above the corresponding convective one while the upper one 
descends significantly below its convective counterpart. As a result, the extension of the absolute 
HMRI beyond the Rayleigh line is considerably shorter than that of the convective instability. The 
absolute HMRI is supposed to be self-sustained and, thus, experimentally observable without any 
external excitation in a system of sufficiently large axial extension. 

PACS numbers: 47.20. Qr, 47.65.-d, 95.30.Lz 



INTRODUCTION 



The magnetorotational instability (MRI) is known to 
be able to destabilize hydrodynamically stable flows by 
means of an externally imposed magnetic fleld as orig- 
inally shown by Velikhov [l| and analyzed in more de- 
tail by Chandrasekhar [Ij for cylindrical Taylor-Couette 
flow of a perfectly conducting fluid subject to an axial 
magnetic field. Three decades later Balbus and Hawley 
01 suggested that, in a similar way, the hydrodynami- 
cally stable Keplerian velocity distribution in accretion 
disks could be rendered turbulent by the MRI account- 
ing for the formation of stars and entire galaxies pro- 
ceeding much faster than it could be accomplished by 
the viscous angular-momentum transport alone. Mean- 
while this proposition has triggered not only numerous 
theoretical and numerical studies 01 but also some ex- 
perimental efforts as well [H, However, one of the 
main technical challenges to laboratory MRI is the mag- 
netic Reynolds number Rm which is required to be ~ 10 
at least. For a liquid metal with the magnetic Prandtl 
number Pm 10~^ — 10~^ this translates into a hydro- 
dynamic Reynolds number Re = Rm/Pm ~ 10^ — 10^ 
Thus, the base flow on which the MRI is supposed to 
be observable may easily be turbulent at such Reynolds 
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numbers independently of MRI as in the experiment of 
Sisan et al. A way to circumvent this problem was 
proposed by Hollerbach and Riidiger Q who suggested 
that MRI can take place in the Taylor-Couette flow at 
Re ~ 10'^ when the imposed magnetic fleld is heHcal 
rather than purely axial as in the classical case. The 
theoretical prediction of this new type of helical MRI 
(HMRI) was soon succeeded by a conflrming experimen- 
tal evidence provided by the so-called PROMISE facility 
Nevertheless, these experimental observations 
have subsequently been questioned by Liu et al. who 
flnd no such instability in their inviscid theoretical anal- 
ysis of flnite length cyHnders with insulating end caps. 
They suspect the observed phenomenon to be a transient 
growth rather than a self-sustained instability (isl. fl^. 

Indeed, such an interpretation of the HMRI is possible 
when the analysis is based only on the conventional linear 
stability analysis for separate Fourier modes as done by 
Hollerbach and Riidiger @| following the classical MRI 
approach. However, there is a principal difference be- 
tween the classical and the helical MRIs, namely, the 
former is stationary whereas the latter is traveling. It is 
important to emphasize that the conventional stability 
analysis for traveling waves yields the so-called convec- 
tive instability threshold at which the system becomes 
able to amplify certain externally excited perturbations. 
At this threshold the perturbation grows in time only 
in the frame of reference moving with its group veloc- 
ity while it asymptotically decays in any other frame of 



reference including the laboratory one [lB|. Eventually, 
such a growing while traveling perturbation reaches the 
end wall where it is absorbed unless the system is able to 
reflect it back. The latter supposes reflection symmetry 
in the system which, however, is not the case provided 
that the magnetic field is helical. Thus, it is indeed un- 
clear whether the HMRI can be self-sustained in an ideal 
Taylor-Couette flow of large but flnite axial extension. 

This question is addressed in the second part of the 
present study where the absolute HMRI is found to exist 
besides the convective one which, in turn, is analyzed in 
detail in the flrst part. Note that the existence of absolute 
instability is nontrivial as known, for instance, for the 
Ponomarenko dynamo which has a convective but no 
absolute instability threshold The latter requires an 
additional return flow to be included in the original Pono- 
marenko model The distinction between convective 
and absolute instabilities is relevant mainly for open flows 
and unbounded geometries [l^. In flnite geometries, it is 
important to distinguish transiently growing and noise- 
sustained perturbations from the self-sustained linear in- 
stabilities, which are always global with the threshold 
asymptotically approaching from above that of the abso- 
lute instability as the system size increases (20l. [2l|. 

We consider both the convective and the absolute HM- 
RIs in the inductionless approximation corresponding to 
Pm = that was suggested in our previous work [22l |. 
This approximation, which leads to a significant simpli- 
fication of the problem, allows us to focus exclusively 
on the HMRI because it does not capture the conven- 
tional MRI We show that the HMRI is effective 
only in a relatively narrow range of the ratio of rotation 
rates of the inner and outer cylinders beyond the limit of 
purely hydrodynamic instability defined by the so-called 
Rayleigh line. For the convective HMRI, the range of 
instability is considerably larger for perfectly conducting 
cylinders than that for insulating ones. In addition we 
find that the HMRI is effective only in a limited range 
of Reynolds numbers. Namely, for any unstable mode, 
there is not only a lower critical Reynolds number by ex- 
ceeding which the HMRI sets in but also an upper one 
by exceeding which it disappears again. It is this upper 
threshold that distinguishes HMRI from a magnetically 
modified Taylor vortex fiow. Absolute HMRI exists in a 
significantly narrower range of parameters than the con- 
vective one. In contrast to the convective HMRI, the 
absolute one is much less dependent on the conductivity 
of the boundaries. 

The paper is organized as follows. In Sec. Ullwe formu- 
late the problem using the inductionless approximation. 
Numerical results concerning the convective and absolute 
instability thresholds for both insulating and perfectly 
conducting cylinders are presented in Sees. IIII Al and 
IIIIB| respectively. Section |IV] concludes the paper with 
a summary and a comparison with experimental results 
of Stefani et al. 
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Figure 1: Sketch to the formulation of the problem. 



II. PROBLEM FORMULATION 

Consider an incompressible fiuid of kinematic viscos- 
ity V and electrical conductivity a filling the gap between 
two infinite concentric cylinders with inner radius Ri and 
outer radius Ro rotating with angular velocities fli and 
flo, respectively, in the presence of an externally imposed 
steady magnetic field Bg — B^e^ -\- B^e^ with axial and 
azimuthal components B^ = Bq and Bff, — PBoRi/r in 
cylindrical coordinates {r,4>,z), where /3 is a dimension- 
less parameter characterizing the geometrical helicity of 
the field (Fig. [l]). Further, we assume the magnetic field 
of the currents induced by the fiuid fiow to be negligible 
relative to the imposed field. This corresponds to the 
so-called inductionless approximation holding for most 
of liquid-metal magnetohydrodynamics characterized by 
small magnetic Reynolds numbers Rm = (iqctvoL <^ 1, 
where /iq is the magnetic permeability of vacuum and 
vq and L are the characteristic velocity and length scale. 
The velocity of fiuid fiow v is governed by the Navier- 
Stokes equation with electromagnetic body force 

-^ + (v-V)v = — Vp + t^V2v+ -j X Bo, (1) 
ot p p 

where the induced current follows from Ohm's law for a 
moving medium 

j = a(E + vxBo). (2) 

In addition, we assume the characteristic time of velocity 
variation to be much longer than the magnetic diffusion 
time To ;2> t„i — poaL^ that leads to the quasi-stationary 
approximation, according to which V x E = and E = 
— V<i>, where $ is the electrostatic potential. Mass and 
charge conservations imply V • v = V • j = 0. 

The problem admits a base state with a purely az- 
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imuthal velocity distribution vo(r) = e0t;o(r), where 
, , ^oRo ~ ^iRf 1 ^^o — 

Note that the magnetic field does not afi^ect the base fiow 
because it gives rise only to the electrostatic potential 
*i'o('') — Bq J Vf){r)dr whose gradient compensates the 
induced electric field, so that there is no current in the 
base state (jg = 0). However, a current may appear in a 
perturbed state, 

{ri}'-'>-{ri:}«+{;;i;}<-" 

where Vi, pi, jj^, and $1 present small-amplitude pertur- 
bations for which Eqs. H]) and ^ after linearization take 
the form 

t^v 1 1 

^ + (vi-V)vo+(vo-V)vi = — Vpi+vV\i + -i^-KBa, 
ot p p 

(3) 

Jl = cr(- V$i + VI X Bo). (4) 

In the following, we focus on axisymmetric perturbations 
which are typically much more unstable than nonaxisym- 
metric ones (Hi- For such perturbations the solenoidity 
constraints are satisfied by meridional stream functions 
for fiuid fiow and electric current as 

V = ve^j, + V X i^l^e^), j = je^ + V x (he^). 

Note that h is the azimuthal component of the induced 
magnetic field which is used subsequently instead of <I> for 
the description of the induced current. Thus, we effec- 
tively retain the azimuthal component of the induction 
equation to describe meridional components of the in- 
duced current while the azimuthal current is explicitly 
related to the radial velocity. The use of the electrostatic 
potential $, which provides an alternative mathematical 
formulation for the induced currents in the inductionless 
approximation, would result in slightly more complicated 
governing equations. In addition, for numerical purposes, 
we introduce also the vorticity uj = we^ + V x (ve^) = 
V XV as an auxiliary variable. The perturbation is sought 
in the normal mode form 

, /ii} (r, t) = {«, ^, h] (r)e-'*+*'=^ , (5) 

where 7 is, in general, a complex growth rate and k is 
the axial wave number which is real for the conventional 
stability analysis and complex for absolute instability. 
Henceforth, we proceed to dimensionless variables by us- 
ing Ri, Ri/vj Ri^i, Bq, and aBoRi^li as the length, time, 
velocity, magnetic field, and current scales, respectively. 
The nondimensionalized governing equations then read 
as 

ji, = DkV + Reik{r'^Qyr^^i]j + }la^ikh, (6) 
7CJ = Dki:j + 2Reiknv-'Ha^ikiikiij + 2l3r-^h), (7) 
= Dk^J + io, (8) 
= Dkh + ik{v -2Pr-^iP), (9) 



where Dkf = r ^ irf')' — [r + k'^)f and the prime 
stands for d/dr, Re = Rfili/v and Ha — RiB^^J a / {pv) 
are Reynolds and Hartmann numbers, respectively, and 

is the dimensionless angular velocity of the base fiow de- 
fined by A — Ro/Ri and p, — flo/^i- The boundary 
conditions for the fiow perturbation on the inner and 
outer cylinders at r = 1 and r = A, respectively, are 
V = ip = tp' — 0. Boundary conditions for h on insulat- 
ing and perfectly conducting cylinders, respectively, are 
/i = and {rhy = at r = 1; A. 

The governing equations (l6l)-([9]) for perturbation am- 
plitudes were discretized using a spectral collocation 
method on a Chebyshev-Lobatto grid with a typical num- 
ber of internal points A'' = 32 — 96. Auxiliary Dirichlet 
boundary conditions for uj were introduced and then nu- 
merically eliminated to satisfy the no-sHp boundary con- 
ditions i^' = 0. The electric stream function h was ex- 
pressed in terms of v and 1/; by solving Eq. ([9]) and then 
substituted in Eqs. ^ and ^ that eventually resulted 
in the 2N x 2N complex matrix eigenvalue problem which 
was solved by the LAPACK's ZGEEV routine. 

III. NUMERICAL RESULTS 

A. Convective instability 

In this section, we consider the so-called convective 
instability threshold supplied by the conventional linear 
stability analysis with real wave numbers k, as done in 
most of previous studies [3, US]- Note that at the 
convective instability threshold the system becomes able 
to amplify certain perturbations which however might be 
not self-sustained and, thus, experimentally unobservable 
without an external excitation. The following results con- 
cern the radii ratio of outer to inner cylinder A = 2 and 
we start with insulating cyHnders which form the side 
walls of the system. 

1. Insulating cylinders 

The critical Reynolds number, wave number, and fre- 
quency are shown in Fig. [2] versus the angular velocity 
ratio p of outer to inner cylinder for Hartmann num- 
ber Ha = 15 and various geometrical helicities /3. For 
/? = corresponding to a purely axial magnetic field, 
the critical Reynolds number tends to infinity as p ap- 
proaches the Rayleigh fine pc = A~^ = 0.25 defined by 
d (r^ri) /dr = 0. Thus, for P — 0, the range of insta- 
bility is limited by the Rayleigh line, i.e., p < pc, as 
in the purely hydrodynamic case. For helical magnetic 
fields defined by /? 7^ 0, the instability extends well be- 
yond the Rayleigh line, as originally found by Hollerbach 
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Figure 2: (a) Critical Reynolds number Rec, (b) wave number kc, and (c) frequency ojc versus ^ at Hartmann number Ha = 15 
and various helicities of the magnetic field for insulating cylinders. 



and Riidiger Note that it is this extension of the in- 
stabiUty beyond its purely hydrodynamic limit, that for 
ideal Taylor-Couette flow is defined by the Rayleigh line, 
which constitutes the essence of the MRI. Comparing the 
stability curves presented in Fig. [2l[a) for fixed Ha to 
those of Hollerbach and Riidiger [Sy , which are presented 
for Pm 7^ and variable Ha yielding minimal RCc, there 
are two differences to note. First, the range of instabil- 
ity is limited by a certain /imax, which depends on the 
helicity f3 and the Hartmann number, as shown in Fig. 
fTsFa). Second, the destabiHzation beyond the Rayleigh 
line is effective only in a limited range of Reynolds num- 
bers bounded by an upper critical value which tends to 
infinity as fi approaches the Rayleigh line from the right. 

The origin of the upper critical Reynolds number, by 
exceeding which the flow becomes linearly stable again, 
is illustrated in Fig. [3] showing (a) the Reynolds number 
and (b) the corresponding frequency of marginally stable 
modes versus their wave number k for /? = 5, Ha = 15, 
and fi = 0.27. As seen, the marginal stability curves 
for ^ beyond the Rayleigh line (/i > 0.25) form closed 
loops which collapse at /i = HmnK- Thus, unstable modes 
exist only within limited ranges of wave and Reynolds 
numbers. Obviously, at sufficiently large Reynolds num- 
bers the flow becomes effectively non-magnetic as inertia 



starts to dominate over the electromagnetic forces sup- 
pressing the HMRI. 

Note that the suppression of HMRI at high Re is re- 
lated to the negligibility of the induced magnetic fleld 
as long as Rm <C 1. In this case, the electric current is 
induced only by the velocity perturbation crossing the 
imposed magnetic fleld. In the conventional MRI con- 
versely to the HMRI, also the induced fleld is relevant, 
which crossed by the base flow induces an additional elec- 
tric current. These two effects are easily noticeable in 
the induction equation. Thus, in the HMRI, the result- 
ing electromagnetic force is not affected by the base flow 
which is not the case for the conventional MRI, where the 
electromagnetic force remains signiflcant with respect to 
inertia also at high Re. 

As seen in Figs. [D^a) andSJb), the critical Reynolds 
number can vary with the Hartmann number in three 
different ways depending on /i. For /i = 0, the critical 
Reynolds number is bounded at Ha = because a purely 
hydrodynamic instability is possible before the Rayleigh 
line. Numerical results evidence that the increase in the 
Hartmann number results in the growth of the critical 
Reynolds number with asymptotics ~ Ha. For fi = 0.25, 
which lies exactly on the Rayleigh line, the flow is hy- 
drodynamically stable without the magnetic fleld. Thus, 
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Figure 3: Neutral stability curves: (a) marginal Reynolds number and (b) the frequency versus the wave number for /3 = 5 and 
Ha = 15 at various ratios /i of angular velocities of the outer to inner cylinder. 



in this case, the critical Reynolds number increases as 
- Ha"^ as Ha ^ because there is no finite value of 
the critical Reynolds number without the magnetic field. 
The corresponding critical wave number tends to a finite 
value independent of (3. With increase in the Hartmann 
number, the critical Reynolds number attains a minimum 
at Ha 10 and starts to grow at larger Hartmann num- 
bers similarly to the previous case. The corresponding 
critical wave number decreases asymptotically as ~ Ha~ 
that means a critical wavelength increasing directly with 
the magnetic field strength. The critical frequency plot- 
ted in Fig. m changes from a constant value of 30 at 
small Ha to another nearly constant value of about 100 
slightly varying with (3 at large Ha. At large helicities 
(/3 = 15), another relatively short-wave instability mode 
dominates up to a Hartmann number Ha ~ 30, where 
the most unstable mode switches back to the long-wave 
one which is characteristic for smaller helicities. Tran- 
sition to this large-/? mode is also obvious in Fig. [2] for 
(3=lbai 0.235. As seen in Fig. B^b) for = 0.27, 
which is beyond the Rayleigh line, there is no instability 
as Ha 0. Consequently, a finite minimal value of Ha 
depending on (3 is necessary in this case. Moreover, the 
instability is limited by the upper branch of the critical 
Reynolds number discussed above which merges with the 
lower branch at the minimum of the Hartmann number 
for the given helicity fi. 

The variation of the critical Reynolds number with the 
heHcity (3 shown in Fig. [Sja) for ~ 0.25 = /ic lying ex- 
actly on the Rayleigh fine differs considerably from the 
other case with /i = 0.27 > /Zc [see Fig. [lb]. In the first 
case, the fiow can be destabilized by the magnetic field 
of however small helicity /? ^ that results in the crit- 
ical Reynolds number increasing as ~ 1//3. For > ^c, 
a certain minimal helicity depending on the Hartmann 
number is needed. Moreover, in this case, there is also 
an upper critical Reynolds number. In both cases, there 
is some optimal /3 « 5 — 8 at which the lower critical 
Reynolds attains a minimum. Further increase in (3 re- 



sults in the growth of the critical Reynolds number with 
a significantly different asymptotic behavior in both con- 
sidered cases. For fi > fic, there is a maximal (3 depend- 
ing on the Hartmann number at which the upper and 
lower branches of the critical Reynolds number merge to- 
gether and the instability disappears whereas there seems 
to be no such merging point at any finite /3 when ^jl = ^c- 
The critical wave number plotted in Figs. Wic) andlUJd) 
is seen to increase with (3 with some jumps at larger Ha 
as discussed above. 



2. Perfectly conducting cylinders 

For perfectly conducting cyHnders, the marginal stabil- 
ity curves shown in Fig. [6] differ considerably from those 
for insulating walls (see Fig. [3]). Although in both cases 
beyond the Rayleigh line large wave numbers (fc ;2> 1) 
are always stable, the range of instability for perfectly 
conducting cyHnders at moderate (3 <\Q extends to ar- 
bitrary small wave numbers fc ^ whereas for insulating 
cylinders it is limited to sufficiently large k. As in the in- 
sulating case, for each unstable mode there is not only the 
lower but also the upper marginal Reynolds number both 
increasing as ~ 1/fc toward small k. Thus, the increase 
in the Reynolds number results in the shift of instability 
to smaller wave numbers, i.e., longer waves. As a result, 
there is no upper critical Reynolds number for moderate 
/3 < 10 when both cylinders are perfectly conducting. 

Consequently, as seen in Figs. (iKa) and[7jb), the crit- 
ical Reynolds number becomes very large while the crit- 
ical wave number tends to zero as ^ approaches some 
critical ^max which varies with (3. The critical frequency 
LUc shown in Fig. [7l[c) tends, respectively, to some fi- 
nite value. This behavior changes at larger (3 becoming 
similar to that for insulating cylinders. As seen in Fig. 
[Tja), for /3 > 10, the curves of the critical Reynolds num- 
ber start to bend back at ^max toward smaller ^ rather 
than tend to infinity. The corresponding critical wave 
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Figure 4: [(a) 
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and (b)] Critical Reynolds number Rcc, [(c) and (d)] wave number kc, and [(e) and (f)] frequency tOc versus 
number for (a) /i = 0, [(a), (c), and (e)] 0.25, and [(b), (d), and (f)] 0.27 at various helicities /3 and insulating 



numbers remain finite whereas the critical frequency in- 
creases with the Reynolds number [see Figs. jTb and|7h]. 
At intermediate /? the Hmiting value of /Zmax, up to which 
the instability extends beyond the Rayleigh line, is seen 
in Fig. flSlfa) to attain a maximum which is consider- 
ably larger than that for insulating walls. At larger /3, 
the limiting values /J,max decrease approaching those for 
insulating walls. 

The dependence of the critical Reynolds number on the 
Hartmann number plotted in Fig. |8l[a) at various helici- 



ties is similar to that for insulating cylinders. First, Rec 
attains a minimum at Ha =7—10 and a finite minimal 
value of the Hartmann number is required for instabil- 
ity when ^ > ^c- For moderate /S < 10, in contrast to 
the insulating case, RCc and kc tend to infinity and zero, 
respectively, as the Hartmann number approaches this 
minimal value which depends on /3. For /3 > 10, the criti- 
cal Rec has an upper branch which merges with the lower 
one at the minimal value of Ha as in the case of insulat- 
ing cylinders. At sufficiently large Hartmann numbers. 
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Figure 5: [(a) and (b)] Critical Reynolds number Rec, [(c) and (d)] wave number fcc, and [(e) and (f)] frequency versus /3 
for [(a), (c), and (e)] /i — 0.25 and [(b), (d), and (f)] /i = 0.27 at various Hartmann numbers. 



the instability is seen to switch to a long-wave mode with 
the critical wave numbers and Reynolds numbers varying 
asymptotically as ^ Ha~^ and ~ Ha, respectively. 



B. Absolute instability 

In this section, we turn to the absolute instability for 
which the wave number k is in general a complex quan- 
tity with real and imaginary parts and fc^, respectively 
[2a, 129(1 . It is important to realize that the convective in- 



stability threshold considered above is not sufficient for 
the development of a self-sustained instability unless the 
system is mirror symmetric along the direction of propa- 
gation which, however, is not the case when the magnetic 
field is helical. The convective instability just ensures 
the ability of the system to amplify external perturba- 
tions excited with the critical frequency. From the math- 
ematical point of view, the problem is that in an axially 
bounded system the perturbation has to meet certain 
boundary conditions at two end walls that, however, can 
not be accompHshed by a single Fourier mode. When 
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Figure 7: (a) Critical Reynolds number Rec, (b) wave number kc, and (c) frequency uJc versus for Hartmann numbers Ha = 15 
and various magnetic field helicities (perfectly conducting cylinders). 



the critical Fourier mode is replaced with a correspond- 
ing wave packet of a limited spatial extension, one finds 
such a perturbation to grow only in the frame of ref- 
erence traveling with its group velocity while it decays 
asymptotically in any other frame of reference including 
the laboratory one which is at rest. The growth of a 
perturbation in the laboratory frame of reference is en- 



sured by the absolute instability threshold, at which the 
group velocity of the wave packet becomes zero. Thus, 
formally, the absolute instability requires one more con- 
dition to be satisfied, i.e., zero group velocity, by means 
of an additional free parameter — the imaginary part of 
the wave number. Note that although the group veloc- 
ity in non-conservative media is, in general, a complex 
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quantity, for the most dangerous perturbations satisfy- 
ing d^r/dkr — 0, it is real and, thus, coincides with its 
common definition. 

Alternatively, the absolute instability may be regarded 
as an asymptotic case of the global instability when 
the axial extension of the system becomes very large 
[isl . m, HBl- The basic idea is that for a convectively 
unstable perturbation to become self-sustained a feed- 
back mechanism is needed which could transfer a part 
of the growing perturbation as it leaves the system back 
to its origin. Such a feedback can be provided by the 
reflections of perturbation from the end walls or, gener- 
ally, by the end regions where the base state becomes 
axially non-uniform. If the base state is both station- 
ary and axially uniform, the coefficients of the linearized 
perturbation equations do not depend on time and on 
the axial coordinate, respectively. Then, as for linear 
differential equations with constant coefficients, the par- 
ticular solution for the perturbation varies exponentially 
in both time and axial coordinate as supposed by Eq. 
([5]) where both the growth rate 7 and the wave number 
k may be in general complex. At the end walls, where 
the base state is no longer axially invariant, the particu- 
lar solutions with different wave numbers become linearly 
coupled while their time variation remains unaffected as 



long as the base state is stationary. Thus, the reflec- 
tion of a perturbation by the end wall in general couples 
modes with different wave numbers but with the same 7. 
Sufficiently far away from the end walls the reflected per- 
turbation is dominated by the mode with the imaginary 
part of the wave number ki corresponding to either the 
largest growth or lowest decay rate along the axis. Con- 
sequently, sufficiently away from the end walls a global 
mode is expected to consist of two such waves coupled 
by reflections from the opposite end walls. Taking into 
account that the amplitude of the reflected wave is pro- 
portional to that of the incident wave, it is easy to flnd 
that in a sufficiently extended system both waves must 
have the same imaginary part ki of the wave number 
whereas the real parts may be different. Additionally, 
for two such waves to be coupled by reflections from the 
end walls, they have to propagate in opposite directions. 



1. Insulating cylinders 

We search for such a pair of modes by considering the 
conventional neutral stability curves at various ki. As 
seen in Figs. [DJa) and (b), the increase in A;,, on one 
hand, results in the reduction in the wave number range 
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Figure 9: Neutral stability curves: (a) marginal Reynolds number and (b) the frequency versus the wave number, and (c) the 
Reynolds number versus the frequency for — 0.26, /3 = 5, and Ha = 15 at various imaginary parts of wave number ki. 



admitting such neutrally stable modes. On the other 
hand, the lower branch of marginal Re and the corre- 
sponding frequency first increase with ki in the whole 
wave number range and then start to decrease at larger 
kr when ki becomes sufficiently large [ki > 1.8). How- 
ever, more important information is obtained by plotting 
the marginal Re and frequency from the previous curves 
against each other as in Fig. [Djc). First, similarly to 
the previous curves, these ones also form closed loops 
that shrink as ki is increased. It is important to notice 
that at sufficiently large ki these loops start to intersect 
themselves in some point as shown in the inset at the top 
of Fig. M,c). This self-intersection, which is of primary 
importance here, occurs only in the limited range of posi- 
tive ki . The point of intersection means that at the given 
Reynolds number there are two modes with the same 
frequency and the same imaginary but possibly different 
real parts of the wave number. As discussed above, two 
such modes could be coupled by refiections from the end 
walls and thus form a neutrally stable global mode in 
an axially bounded system provided that they propagate 
in opposite directions (25l . \2^ . To determine the direc- 
tion of propagation we use a local criterion [23| which 
we showed to be equivalent to the Briggs pinching crite- 
rion [Hi for the upper instability branch at the given ki. 



Namely, the direction of propagation of both intersecting 
branches can be deduced from their variation with ki. If 
upon a small variation of ki one branch rises to higher Re 
while the other descends to lower Re, which is the case 
here, it can be shown that both intersecting branches 
correspond to oppositely propagating modes. The low- 
est possible Reynolds number admitting two such modes 
is attained when the loop below the intersection point 
collapses to a cusp as seen in the inset at the bottom of 
Fig. [9^c) for ki = 1.8. The cusp is formed as both in- 
tersection points of the loop merge together. It means 
that at the cusp point not only the imaginary but also 
the real parts of both wave numbers become equal. This 
point corresponds to the absolute instability at which the 
length of the wave packet of the global mode formed by 
two waves with merging wave numbers tends to infin- 
ity. Further we focus on this absolute instability which, 
in contrast to the convective one considered above, can 
be self-sustained in a sufficiently extended system. Note 
that the approach outlined above to find the absolute in- 
stability is an extension of the well-known cusp r nap for 
the complex frequency plane to the (Re-o;) plane (29l . [30| 
by using the neutral stability condition Xr (Re) = which 
maps the real part of the growth rate to the marginal 
Reynolds number. 
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Figure 10: (a) Critical Reynolds number Rec, (b) frequency uic, (c) real and (d) imaginary parts of wave number kr,c and ki, 
versus ^ at Ha = 15 and various magnetic field helicities for insulating cylinders. 



The critical Reynolds number, frequency, and the crit- 
ical complex wave number for the absolute instability 
threshold is plotted in Fig. [lO] versus ^ at Ha = 15 
and various helicities of the magnetic field. Comparison 
virith the corresponding convective instability, the crit- 
ical parameters of which are plotted in Fig. [21 shows 
that before the Rayleigh Hne the threshold of absolute 
instability is only sHghtly above the convective one. The 
difference between both thresholds becomes significant at 
the Rayleigh line. Although the absolute instability ex- 
tends beyond the Rayleigh Hne when the magnetic field is 
heHcal, the range of extension is noticeably shorter than 
that of the convective instability (see Fig. [13]). More- 
over, the upper critical Reynolds number for the absolute 
instability is considerably lower than that of the convec- 
tive one. Although the difference between the critical 
Reynolds numbers for the absolute and convective insta- 
bilities is insignificant before the Rayleigh line, the crit- 
ical wave numbers for the absolute instability shown in 
Fig. \Wic) are considerably larger than those for the con- 
vective instability [see Fig. [2b]- This difference increases 
with P that results in the rise of the critical wave number 
for absolute instability while the increase in the corre- 
sponding quantity for the convective instability thresh- 
old is insignificant. In contrast to this, the imaginary 



part of the critical wave number for the lower instabil- 
ity branch fc^ « 1.8 is almost invariable with both /? and 
/i except for (3=1, where a jump of the instability to 
a larger wave number takes place at w 0.235. Note 
that positive ki corresponds to the amplitude of the crit- 
ical perturbation growing axially downward which is an 
additional feature predicted by the absolute instability. 
Beyond the Rayleigh line the absolute instability simi- 
larly to the convective one is effective only in a limited 
range of Reynolds numbers which is bounded from above 
by the upper critical branch tending to the Rayleigh Hne 
from the right as the Reynolds number increases. The 
critical complex wave number is seen in Figs. [TOlfc) and 
[TOF d') to tend to a certain limiting value independent of 
/?. 



2. Perfectly conducting cylinders 

The neutral stability curves for perfectly conducting 
cylinders plotted in Fig. [11] are seen to start forming 
closed loops when fc^ > becoming similar to the corre- 
sponding curves for insulating cyHnders shown in Fig. 
[9] In a certain range of ki the curves of the critical 
Reynolds number plotted against the frequency in Fig. 
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Figure 11: Neutral stability curves: (a) marginal Reynolds number, (b) the frequency versus the wave number, and (c) the 
Reynolds number versus the frequency for ^ — 0.26, /3 = 5 and Ha — 15 at various imaginary parts of wave number ki for 
perfectly conducting cylinders. 



fUT c) intersect themselves that imphes the existence of 
two neutrally stable modes with the same Reynolds num- 
ber, frequency, and imaginary part of the wave number 
but different real parts of the wave number. As discussed 
above, two such modes can be coupled by reflections from 
the end walls and thus form a neutrally stable small- 
amplitude global mode in the system of a large but finite 
axial extension provided that those modes propagate in 
opposite directions that is implied by the variation of the 
marginal Reynolds number upon a small variation of ki . 

For perfectly conducting cylinders, the critical 
Reynolds number, frequency, and complex wave number 
for the absolute instability threshold plotted in Fig. [12] 
versus ^ for Ha = 15 and various helicities differ signifi- 
cantly from the corresponding critical parameters for the 
convective instability threshold (see Fig. [7|). First, the 
range of extension of the absolute instability beyond the 
Rayleigh line is much shorter than that of the convec- 
tive instability. Note that in contrast to the convective 
instability there is no significant difference with respect 
to the extension of the absolute instability beyond the 
Rayleigh line between insulating and perfectly conduct- 
ing cylinders (see Fig. [13]). Second, beyond the Rayleigh 
line, similarly to insulating cylinders, for all (3 the range 



of unstable Reynolds numbers is bounded from above by 
the upper critical branches which approach the Rayleigh 
line from the right as the upper critical Reynolds number 
tends to infinity. Similarly to the insulating cylinders, the 
corresponding critical complex wave number tends to a 
certain asymptotic value independent of (3 [see Figs. [T2b 
and [T2ti] . Third, beyond the Rayleigh line the critical 
wave numbers for the absolute instability are noticeably 
greater than those for the convective instability, espe- 
cially for /? ^ 10 when the critical wave numbers for the 
convective instability tend to zero [see Fig. Wp\- 



IV. CONCLUSION 

In this study we have analyzed numerically the MRI 
of Taylor-Couette flow with a helical external magnetic 
field. The problem was considered in the inductionless 
approximation defined by a zero magnetic Prandtl num- 
ber (Pm — 0). First, we carried out a conventional linear 
stability analysis for perturbations in the form of Fourier 
modes specified by real wave numbers. The helical mag- 
netic field was found to extend the original instability 
to a relatively narrow range beyond its purely hydrody- 
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Figure 12: (a) Critical Reynolds number Rec, (b) frequency ojc, (c) real and (d) imaginary parts of wave number kr,c and ki^c 
versus ^ for Ha = 15 and various magnetic field helicities /3 for perfectly conducting cylinders. 




namic limit defined by the Rayleigh line. The range of 
destabilization was found to be considerably larger for 
perfectly conducting cylinders than that for insulating 
ones. For insulating cylinders, the instability beyond the 
Rayleigh line is effective only in a limited range of wave 
and Reynolds numbers. Unstable Reynolds numbers are 
bounded by an upper critical value which tends to in- 



finity right beyond the Rayleigh line. For perfectly con- 
ducting cylinders and moderate helicities of the magnetic 
field, the range of unstable wave numbers is bounded only 
from the short-wave end. Although there is an upper 
marginal Reynolds number for each unstable wave num- 
ber, no bounded upper critical Reynolds number exists 
in this case because the range of unstable wave numbers 
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extends to zero, i.e., infinitely long waves. Nevertheless, 
at sufficiently large helicities, the range of unstable wave 
numbers becomes bounded also from below, and an up- 
per critical Reynolds number appears in the same way as 
for insulating cylinders. 

It is important to note that these instabilities predicted 
by the conventional stability analysis in the form of single 
traveling waves correspond to the so-called convective in- 
stability threshold at which the system becomes able to 
amplify certain externally imposed perturbations that, 
however, are not self-sustained and thus may be exper- 
imentally unobservable without a proper external exci- 
tation. The problem is that convectively unstable per- 
turbations grow asymptotically in time only in the frame 
of reference traveling with their group velocity, whereas 
they decay in any other frame of reference including the 
laboratory one. For an instability to be self-sustained and 
thus observable it has to grow in the laboratory frame of 
reference. In an extended system, this condition is sat- 
isfied by the so-called absolute instability which ensures 
a zero group velocity of a growing perturbation. This 
additional condition is satisfied by regarding the wave 
number as a complex quantity with a nonzero imaginary 
part which describes an exponential axial modulation of 
the wave amplitude. Using this concept, we found that 
there is not only a convective but also an absolute HMRI 
implying that this instability can be experimentally ob- 
servable in a system of sufficiently large but finite axial 
extension. In the hydrodynamically unstable range be- 
fore the Rayleigh line, the threshold of absolute insta- 
bility is slightly higher than the convective one. Nev- 
ertheless, the critical wavelength for absolute instability 
is significantly shorter than that for the convective one 
that may allow us to distinguish between both. The ab- 
solute instability threshold rises significantly above the 
convective one beyond the Rayleigh line. As a result, the 
extension of the absolute instability beyond the Rayleigh 
line is considerably shorter than that of the convective 
instability without a marked difference between insulat- 



ing and perfectly conducting cylinders in contrast to the 
convective HMRI. 

The extension of HMRI beyond the Rayleigh line is of 
particular interest from the astrophysical point of view 
regarding a Keplerian velocity profile |l2l . l31| . For a 
Couette- Taylor fiow with a radius ratio A = 2 considered 
here, the Keplerian velocity profile approximately corre- 
sponds to a ratio of rotation rates of — A"^/^ ~ 0.35. 
As seen in Figs. [TSlfb) and [14] for the absolute insta- 
bility, no such value of ^ is reached up to /? = 30 and 
Ha = 150. Whether or not it can be reached at higher f3 
and Ha is still an open question requiring a more detailed 
study using either higher numerical resolution or asymp- 
totic analysis. On the other hand, HMRI in a system of 
large axial extension with a radius ratio of A = 2 might 
be of limited astrophysical relevance for accretion disks 
anyway. 

Conversely to Liu et al. [i3,[l3| we find that the HMRI 
can be a self-sustained instability rather than just a tran- 
sient growth. This contradiction may be due to a couple 
of additional simplifications underlying the analysis of 
Liu et al. First, the viscosity is neglected. Second, the 
electromagnetic force is treated as a small perturbation 
which is a sensible approach within the inviscid approxi- 
mation where an infinitesimal magnetic field can cause 
a correspondingly slow growth of an unstable Fourier 
mode. However, such a perturbative approach may be 
inadequate for the absolute instability which requires a 
finite temporal growth rate of the corresponding convec- 
tively unstable Fourier mode. In addition, note that our 
results do not support the recent findings of Liu [33 |. 
According to his estimates for the cases considered here, 
the absolute HMRI occurs above Re ~ 10^ only. We pre- 
dict the absolute HMRI to occur at the lower value of 
Re 10^, and to disappear again above Re ~ 10^ due 
to the mechanism discussed in Sec. IIII A 11 This dis- 
agreement may be due to the absolute instability analy- 
sis which is carried out by Liu using real wave numbers 
only. Accor ding to the conventional absolute instability 
analysis jH, [la], such an approach yields the long-time 
asymptotics only in the frame of reference traveling with 
the group velocity of the fastest growing perturbation. 
Moreover, owing to its linearity, the analysis is limited to 
sufficiently small perturbation amplitudes only. Namely, 
it is limited up to the point where the first exponentially 
growing perturbation appears, which is the convective 
instability threshold. 

Finally, let us compare the convective and absolute 
instability thresholds calculated for perfectly conducting 
cylinders with the experimental data of Stefani et al. [ll| 
who reports on the observation of HMRI-like traveling 
waves at ^ = 0.27, Re =1775, and a fixed rod current 
6 kA for the coil currents 40 — 100 A that corresponds 
to (3 K, 7.4-3 and Ha w 6.3-15.8. Another observation 
was done at the same /j, but different other parameters: 
Re = 1479, /3 = 6, and Ha = 9.5. As seen in Fig. [HI 
in both cases the experimental points lie well inside the 
range of ^ for convective instability but outside that for 



15 



3000 

^„ 2500 

S3 
n 

1 2000 
« 

2 1500 
o 

c 

^ 1000 

a 
o 

S 500 

S-H 

J 





1 — ■ r 

convective — 
absolute 



0.05 




P = 3, Ha= 15.8 
P = 4.5, Ha= 12 
P = 7.4, Ha = 6.3 



Pi 



3000 
2500 
2000 
1500 
1000 
500 




1 1 r 

convective 

absolute 



P = 6, Ha = 9.5 



0.05 



0.1 




0.15 



0.2 0.25 0.3 



Figure 15: Comparison of the critical Reynolds number Rec for convective and absolute instability thresholds in the case of 
perfectly conducting cylinders with the experimental values (dots) for which traveling waves have been observed at (a) /i = 0.27, 
Re = 1775, /3 « 7.4 - 3, and Ha « 6.3 - 15.8 and (b) Re = 1479, /3 = 6, and Ha = 9.5. 



absolute instability. This discrepancy with the experi- 
mental observations may be due to the deviation of the 
real base flow from the idealized one used in this study. In 
particular, the Ekman pumping driven by the end walls 
in the experiment, which is not taken into account in the 
present analysis, may affect the hydrodynamic stability 
limit of the base flow, i.e., its actual Rayleigh line, which 
however serves as the reference point for the observation 
of MRI. A more detailed comparison with the experi- 
mental observations lies outside the scope of the present 
paper. 

In conclusion, the main result of the present paper is 
the finding of absolute HMRI in addition to the con- 
vective one which can be self-sustained and, thus, ex- 



perimentally observable without external excitation in a 
system of sufficiently large axial extension. A charac- 
teristic feature of HMRI is the upper critical threshold 
existing besides the lower one that distinguishes it from 
a magnetically-modified Taylor vortex fiow. 
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